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Abstract 

We investigate potential quantum nonlinear corrections to Dirac's 
equation through its sub-leading effect on neutrino oscillation proba- 
bilities. Working in the plane- wave approximation and in the fj, — r 
sector, we explore various classes of nonlinearities, with or without 
an accompanying Lorentz violation. The parameters in our models 
are first delimited by current experimental data before they are used 
to estimate corrections to oscillation probabilities. We find that only 
a small subset of the considered nonlinearities have the potential to 
be relevant at higher energies and thus possibly detectable in future 
experiments. A falsifiable prediction of our models is an energy depen- 
dent effective mass-squared, generically involving fractional powers of 
the energy. 

1 Introduction 

It appears that constant neutrino masses, though still not directly confirmed, 
are the simplest way of explaining of current data on neutrino oscillations 
p. Other possibilities, such as Lorentz violating dispersion relations |2J, do 
not seem to be possible explanations of the leading order effects [3]. 

Neutrinos are a valuable probe of new physics because of their weak in- 
teractions and in this paper we will study how nonlinear modifications to 
the quantum mechanical propagation of a neutrino, that is a nonlinear Dirac 
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equation, affects neutrino oscillations. The propagating neutrino wavefunc- 
tions are assumed to follow 



{i^d^ -m + F)ij = 0, 



(1) 



where F is a matrix in spinor space, depending on the wavefunction ip, its 
adjoint and their derivatives [I]. One may think of F as probing deviations 
from exact quantum-linearity, a possibility that had had some empirical tests 
in the non-relativistic regime placing bounds on the nonlinearity scale [3 16] . 

In [7] it was suggested that quantum nonlinearities might be related to 
Lorentz violation, leading one to consider higher-energy tests. Although in 
this paper we initially adopt the more general possibility that F might be 
nonlinear but Lorentz invariant, we find that Lorentz invariant nonlinearities 
are not likely to be probed by neutrino oscillations. Furthermore, we leave 
open the possibility that F might only be an effective nonlinearity, summaris- 
ing unknown microscopic physics, rather than a fundamental modification of 
quantum theory. Thus while in our construction of nonlinear Dirac equations 
in [1] we required the equation to be invariant under the scaling ip —>■ Xi/j, as 
in the linear theory [5], thus leading to nonpolynomial nonlinearities, here 
we consider also simpler polynomial nonlinearities which do not have that 
invariance. 

We restrict our study to the \l — t sector as this oscillation is more likely 
to be probed at higher energies in the near future compared to the e — ji 
oscillation [8] . As in the standard formalism we take the weak eigenstates of 
the neutrinos in the fx — r sector to be superpositions of the mass eigenstates, 



where ip a (x) are the neutrino flavor eigenf unctions, ipi(x) are the "oscillating" 
eigenfunctions, and U is the Leptonic mixing matrix. For two-neutrino flavor 
oscillations the mixing matrix is 



The mass eigenstates are assumed to be plane wave solutions of the nonlinear 
Dirac equation flTJ with modified dispersion relations, 



where Uj are spinors. Thus after travelling a distance L between source and 
detector, and averaging over the unobserved travel time, as in [9] , one obtains 
the usual flavour change probability formula 




(2) 






(4) 





where E is the beam energy, the momenta are taken in the direction of L 
and Ap = Pi — Pj. Recall that in the F = case, with small masses and 
large energies, one has Ap = Am 2 /2E where Am 2 = m 2 — m 2 . Maximum 
oscillations occur when L = Lq, where 

L Ap 7i 7rh 

^r = 2 ~* Lo = Ap~- 



(6) 



In the last step, we have restored the ft's and c's. L is the oscillation length, 
the path-length needed for a neutrino of flavour a to maximally oscillate 
to a neutrino of flavour [3. Thus the oscillation length in the conventional 
approach is given by 

Am^c 4 

The expressions (J5J) and (jHJ) are valid even in the nonlinear theory but 
with a modified Ap. Since we are adopting the plane wave approximation, the 
only difference from the conventional formalism will be modified dispersion 
relations in the expression for Ap. 

We classify the types of F in (|T]) into polynomial (P) or the non-polynomial 
(NP) forms studied in [3], then into Lorentz violating (LV) or Lorentz invari- 
ant (LI). In [3] we had studied the NP type of F's in a double expansion in 
the degree of nonlinearity n and number of derivatives d, for example 



(8) 



has d = and degree n. In this paper, for simplicity we consider nonlineari- 
ties F which consist of a single term with n = a and structure F = (*) Q with 
(*) containing at most one derivative, so that F has at most a derivatives. 
Furthermore we consider the general case of real a, not necessarily an inte- 
ger. While non-integral powers might not be surprising in an effective theory, 
interestingly the specific fraction a = 1/3 appears when one demands con- 
formal invariance of a simple nonlinear Dirac equation [TU]. Also, as noted in 
[3], in the nonpolynomial case one can still preserve separability for general 
a through an appropriate construction. 

Since F is a matrix in spinor space, we will consider two special cases, F oc 
/ and F oc 7 M , representing corrections to the mass or kinetic terms of the 
usual Dirac equation. As the nonlinearity must be small on phenomenological 
grounds, we can compute the modified dispersion relations in perturbation 
theory using the plane wave solutions of the linear theory, see [3]. 

We describe the Lorentz violating cases using background gauge fields 
[TT] . A^, but only consider time-like fields = (A , 0) so that there are no 
direction dependent terms in our oscillation probabilities. The background 
fields will simultaneously control the magnitude of both the nonlinearity 
and the Lorentz violation. A more detailed discussion of Lorentz violating 
nonlinear Dirac equations is in [3]. 



The rest of the paper is structured as follows: In the next section we 
consider one example of a nonlinearity from the class NP-LV in detail and list 
results for other cases. In Section (3) we bound our parameters using current 
experimental data and then use our expressions for the modified oscillation 
probabilities to estimate corrections in future higher energy experiments. In 
Section (4) we discuss the other classes of nonlinearities while Section (5) 
concludes the paper. 

The conventions we use are similar to those of the textbook [12j and 
Ref. [I]. We work in 3 + 1 dimensional flat spacetime with metric g^ u = 
(1,-1,-1,-1), set h = 1 = c and restore them as and when needed. 



2 The NP-LV class 



2.1 Fi 

An example of an F of the type considered in j3] is given by 

where is a real constant background field and a is any real number. Thus 
this is a NP-LV type of nonlinearity of degree n — a, with F proportional to 
/ and no derivatives, d = 0. 

Perturbing around the plane wave limit gives F 2 = (^) • The modified 
dispersion relation is then given by 

k 2 = m 2 -2m(—) +0(A 2a ). (10) 
\ m J 

We take the background field to be time-like. This is to avoid the directional 
dependence coming from the spatial part of the background field. Thus we 
have 

'AEy 

and for small masses 



E 2 ~p 2 = m 2 -2m( (11) 

1 m ' 



m 2 . tn (E\ a - x 



p - E -2E +Aa U • ^ 



The momentum difference is given by 

Am 2 „„ , . ( A a 



= -^r - E-'A — (13) 



2E ~ \m a ~ l 

A' 



where A ( ^sbr ) = l_i s=t. Note that we have indicated a possible 

species dependence in the background gauge field. The oscillation length is 

2nE 



Am 2 - 2E a A ( 



(14) 



and may be written in the form 



2-kE 

~ Am 2 (1 - X x ) 

2.2 Summary of Other NP-LV Cases 

We list here the other types of F we study in this class. 



F oc I, n = d — a. 



iA, 



F oc 7 M , n = a, d = 



F 3 = A^Y 



B u 



F oc 7 M , n = d = a 

F 4 = A^Y 



(15) 



(16) 



(17) 



(18) 



All background fields are real, with only a nonzero time component. The 
dispersion relations are obtained by perturbing around plane waves as in the 
previous subsection to obtain the corresponding X's in formula (ITS)) . 



3 Empirical Bounds and Estimates 

We have the modified oscillation lengths in the form 

2tiE 

~ Am 2 (l -X) 



(19) 



where X is the leading order correction depending on the nonlinearity pa- 
rameters and energy. Current neutrino oscillation data fit the conventional 
neutrino mass scenario quite well. Still, as there are the usual experimen- 
tal uncertainties, we use those to estimate the size of X. From an eyeball 
inspection of the data, we conservatively take X to be in the range 10% to 
0.01% and will use this to constrain the range of values which a can take for 
each type of nonlinearity considered. Following that, we will estimate the 
value of X in future higher energy experiments. 



Since the background fields play the dual role of Lorentz violating and 
nonlinearity parameters, we rewrite them as follows 



(20) 
(21) 
(22) 
(23) 

Upon restoring the factors of ft's and c's, the length dimensions of the pa- 
rameters is given by 



For Fi. 


A a - 


-> ei 


For F 2 : 


A a - 




For F 3 : 


AB a - 




For F 4 : 


AB a - 


-> e 4 



L- 1 

L a-1 



(24) 
(25) 
(26) 
(27) 



3.1 Current Experiments 

From current data [13J, we may take the size of the small dimensionless LV 
parameters to be 



./ 



10 



-27 



(28) 



and if the parameters are neutrino species dependent then we assume Af 
f. The relevant neutrino data are taken to be 



Am 2 
E 



2.5 x 10- 3 eV 2 
WOGeV , 



(29) 
(30) 



that is we use the mean beam energy E = lOOGeV. We assume that the 
mass of neutrinos is of the same order of magnitude as V Am 2 . In the same, 
"order of magnitude" spirit, we estimate the expression A (me) by (Am) (e). 
The nonlinear parameter itself may be written, taking F 2 as an example, 



62 



AT 1 / 



(31) 



where A2 is the nonlinearity length scale. As linear quantum mechanics works 
well, we will require the nonlinearity length scale to be much smaller than 
the Compton wavelength of the neutrino, 



Ac 



h 

mc 



(32) 



but larger than the Planck length X p so that we may continue to work in 
the flat space approximation (however note that the quantum nonlinearity 



might plausibly already include, in an effective way, the effects of gravity). 
Actually, as we shall see below, in conjunction with 10~ 4 <I<C 10~ 4 often 
only one of the constraints A <C A c or Xp <C A is actually useful, the other 
being automatically satisfied. 

As the actual mass of neutrinos may be 1 or 2 orders of magnitude larger 
than what we have assumed above, we compensate for this possibility and 
the other approximations above by taking 

X of order KT 1 to 10~ 4 (33) 
of order KT 1 to 10" 4 (34) 



A c 
A_ 

X p 



of order 10 1 to 10 4 (35) 



The various expressions of X's upon restoring the fc's and c's, and using 
the above assumptions, are given by 

*i = ^L,^ (36) 

A 2 (Am 2 c 4 ) (a+1)/2 V ' 

X2 = (tJr^AmZf^ (37) 
X3 = A 4 (AmV/ a+2,/2 m 



2 a + l A a-l A/ ^ 

(hcf- 1 Am 2 c 



X, = 71 , a -L E a+l (39) 



Our procedure is as follows: Since X depends on the characteristic length, 
A, we invert the relationship to plot A as a function of a using the values for 
/, m, E mentioned above and for the chosen range of X values's. From these 
plots, we can determine the values of a's for which the nonlinearity scale A 
lies within the required range. The plots of A versus a for X\ and X 2 are 
shown in Figure 1 and Figure 2 respectively. From those figures, and similar 
ones for the other X's, we obtain the corresponding range of values for the 
en's: 

For Fi 1.6 < en < 2.1 (40) 
For F 2 1.8 <a 2 < 2.4 or 0.2 < a 2 < 0.4 (41) 
For F 3 0.6 <a 3 < 1.1 (42) 
For F A 0.9 <a 4 < 1.1 (43) 

Since we expect the nonlinear effects to be small, the en's are more likely to 
be observed near the lower bound in X. These en's are indicated in boldface 
in the above equations. Note that the bounds depend on the choice of X and 
A c , A p and so must be updated as more accurate data becomes avialable. 



3.2 Future Experiments 



We now reverse the argument. We plot X versus E for the range of a 
values determined in the previous section and for the allowed regions of A as 
indicated in Eqns. fl34ll35l) . From these plots, we can estimate the values of 
X's expected in future experiments where higher energies will be available 

The plots are shown in Figures 3 to 5 for some cases and parameter values. 
The general trend of X(E) can be seen from equations (1361 - 139]) . Of course 
since we evaluated the X's perturbatively, the expressions and plots are valid 
only for small X, while for larger X the indicated trend is only qualitative. 

4 Other Classes of Nonlinearities 
4.1 P-LV 

Polynomial type of nonlinearities lead to a non-separable Dirac equation. For 
example from the NP-LV, F oc / cases considered before one may remove the 
denominators to get corresponding P-LV type of nonlinearities. But now one 
sees that the modified dispersion relation will depend on the normalisation 
chosen for the wavefunctions. If one chooses the usual plane wave normalisa- 
tion such that ipif} = 1 then the results for the P-LV cases mentioned above 
would be the same as for the NP class. 

So to obtain new results let us explore the popular normalisation i/j^i/j = 1 
which makes ipip = m/E. But this energy factor from the normalisation will 
cancel that from the nonlinearity in the two F oc I P-LV cases obtained in 
the previous paragraph, thus making the modified dispersion relation energy 
independent. Therefore only the F oc 7^ cases are relevant. We label these 
as 

F 5 = A^iWY (44) 
F 6 = Arf[iB v ($dr$ (45) 

Note that while F$ is just the numerator of _F 4 , F 5 is not the numerator of F 3 
as now we have a simpler way to generate a P-LV n = a, d = nonlinearity. 
The corresponding X's are 

2hc(Am 2 c 4 ) {a ~ 2)/2 Af 

X 5 = C[ m C ) — E 1 " (46) 

A5 

2"+iAr 1 (AmV) (Q - 2)/2 Af 

6 /f- \a—l \ 1 

(he) 

and the range of a values are found to be 

< a 5 < 0.9 (48) 
0.5 <a 6 < 0.6 (49) 



Again the a's in boldface are near the more likely values. 

4.2 NP-LI 

If the nonlinearity is Lorentz invariant, the modified dispersion relations 
remain covariant E 2 = p 2 + M 2 but with an effective mass M that depends on 
the nonlinearity parameters [I]. If we take these nonlinearity parameters to 
be non-universal, meaning that the different neutrinos get different effective 
mass corrections, then the oscillation probabilities are modified. However as 
the NP type of nonlinearities are invariant to the normalisation of if}, the 
modification is energy independent and thus not relevant for high-energy 
tests. 

4.3 P-LI 

Choosing the same energy dependent normalisation as in the P-LV cases 
discussed above, and noting the discussion in the previous subsection, one 
sees that in P-LI cases such as {ipip) a the X decreases with increasing energy 
and thus becomes irrelevant at high energies. (We have looked at a few 
simple cases of F in the class P-LI and they show a similar trend). 

5 Discussion 

The way we have parametrized our corrections, the modification to the con- 
ventional neutrino oscillation probabilities may be described in terms of an 
effective energy dependent mass-squared difference, 

Am 2 eff (E) = (1 - X)Am 2 (50) 

Hopefully such an effect may be detectable with higher statistics and energies 
in future experiments. Note that our X's are positive because we took the 
constant in front of the F's in (TjQ) to be positive and absorbed it into the 
nonlinearity parameter. More generally then, the right-hand-side of ( 15 Op 
should read (l±X)Am 2 so that the effective mass might increase or decrease 
with energy. 

Among the various types of nonlinearities we have studied, we find that 
only six have the potential to be detected in future higher energy experiments 
through their increasing energy dependent effect on the neutrino oscillation 
probabilities. In the two F oc I cases, F2, each of the discrete symmetries 
is preserved while in the remaining F oc 7* 1 cases the nonlinearities are C 
and CPT odd. Thus the discrete symmetries might be one way of partly 
discriminating among the possibilities. The mildest energy dependence is 
seen for F$, Fq but note that in those cases the nonlinearity is dependent on 
the energy dependent normalisation. 



Since we worked in the plane wave approximation, the above-mentioned 
effects only probe modified dispersion relations rather than the nonlinearity 
itself. However there are various ways of distinguishing our results from 
other proposals in the literature. Firstly, we found that in modeling the 
nonlinearity by a single term F ~ (*) a in the evolution equation, a turns 
out to have generically noninteger values so that 

Am 2 (E) - Am 2 (0) fl , x 

V o/ s — oc ±EP 51 

Am 2 (0) V ; 

for some positive and typically fractional (3. While even this can be ob- 
tained simply from a modified dispersion relation, [Til [15] . independent of 
a quantum nonlinearity, the availability of a quantum evolution equation in 
our approach will enable a more refined analysis of the phenomena when 
sufficient data become available. We also emphasize that our nonlinearity 
simultaneously violates Lorentz invariance and that is an additional distin- 
guishing feature. 

Further work in this direction would involve going beyond the plane wave 
approximation, leading to genuine nonlinear effects and perhaps leading to 
an understanding of the mixing angles too [TJ. Also, a subleading direc- 
tional dependence of the oscillation probability can be examined by choosing 
space-like background fields. Finally, one should explore if current oscillation 
data can be fit using purely energy dependent effective neutrino masses as 
suggested, for example, in [TJ. 
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6 Figures 



Figure 1: This is a plot of A vs a for X\. The vertical axis, plotted in log 
scale, has units of metre while the horizontal axis is dimensionless. The solid 
and the dashed lines represent X = 1CT 1 and X = 1CT 4 respectively. The 
horizontal lines are the bounds 10 _1 A C and 10~ 4 A C . 

Figure 2: This is a plot of A vs a for X 2 . The vertical axis, plotted in log 
scale, has units of metre while the horizontal axis is dimensionless. The full 
and dashed lines represent X = 10 _1 and X = 10~ 4 respectively. The upper 
two horizontal lines are the bounds 10 _1 A C and 10 _4 A C while the lower two 
horizontal lines are the bounds 10A P and 10 4 A P . 

Figure 3: This is the log-log plot of Xi vs energy. The full and dashed 
lines have a values of 1.6 and 2.1 respectively. Here we have assumed Ai = 
0.1 A c . Note that X 3 (E) has an identical plot to Xi(E) after the following 
redefinition of a: The full and dashed lines have a values of 0.6 and 1.1 
respectively. Again we assume A3 = 0.1 A c . 

Figure 4: This is the log-log plot of X2 vs energy. The full and dashed lines 
have a values of 1.8 and 2.4 respectively, where we have assumed A 2 = 10~ 4 A C . 
The dot line has a = 0.2 where we have assumed A 2 = 10A p . Note that for 
a = 0.4, the energy required for X 2 to be significant (say X = 0.2) is much 
higher than the rest(~ lO 11 ^^). Thus it is not plotted. 

Figure 5: This is the log-log plot of X4 vs energy. The full and dashed lines 
have a values of 0.9 and 1.1 respectively. Here we have assumed A 4 = 10 4 A P . 
Note that if we choose A4 = 10 _4 A C , the results will still be of a similar order 
of magnitude: The reason is that the exponent for A 4 is close to zero. 
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